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I . INTRODUCTION; Convex-Profile Inversion 

The initial six months of research sponsored by NASA Grant 
NAGW-543 has proceeded rapidly and has yielded a variety of important 
results. The focus of this research is a new approach, called 
"convex-profile inversion," to the interpretation of asteroid 
lightcurves. This technique, developed by the principal investigator 
in collaboration with R. Connelly, was introduced via the published 
literature this spring and is already being recognized as the most 
important contribution since that of H. N. Russell (1906) to solving 
the problem of extracting information about an asteroid's shape from 

I 

its lightcurve. Except for the relative handful of objects observable 
with radar or during stellar occultations, lightcurves provide the 
only clues to an asteroid's shape. Unfortunately, the form of a 
lightcurve is influenced by the vi ewing/i 11 umi nation geometry and the 
asteroid's light -scattering behavior, as well as by its shape. Even 
if we had lightcurves for an asteroid for all physically possible 
Sun-Earth-asteroid configurations, it would still be impossible to 
determine the asteroid's shape from disk-integrated measurements. All 
we can hope to do is to obtain constraints on the shapes of asteroids. 

Prior to convex-profile inversion (CPI), almost all efforts to 
obtain such constraints have invoked an axi symmetric shape described 
by only one or two parameters. These analyses necessarily discard 
much of the interesting information that distinguishes individual 
lightcurves, including all odd harmonics. 

CPI obtains a convex profile _P from an asteroid's lightcurve. 

The number of parameters that characterize the profile is limited only 
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by the number of Fourier harmonics used to represent the parent 
lightcurve, so our inversion method can preserve much (if not all) of 
the salient information contained in the parent lightcurve. Whenever 
four ideal conditions are Satisfied, P_ is an estimator for the 
asteroid's "mean cross section" £, a convex set defined as the average 
of all cross sections £(z) cut by planes a distance z above the 
asteroid's equatorial plane. £ is therefore a 2-D average of the 
asteroid's 3-D shape. The ideal conditions are: (A) each curve £(z) 

is convex, (B) the asteroid's scattering law is uniform and geometric, 
(C) the astrocentric declinations 6$» 5 E of the Sun and Earth are 
zero and (D) the solar phase angle 4> is known and nonzero. (CPI's 
geometrical conventions are such that for any given Sun-Earth-asteroid 
configuration, the sign of the solar phase angle will correspond to 
either retrograde or direct rotation.) 

To obtain a profile from a lightcurve, we first find the 
lightcurve's Fourier series: 1(9) = l Cpe 1119 and define the 

radius of curvature function: r(e) = j c^e 1 ’ 09 , where 

dji = Cji/vjn (1) 

and _Vn is a known function (Fig. 1) of the harmonic index n and the 

solar phase angle <j>. The profile's cartesian coordinates are 

( -/® r(t) sin t dt, /® r(t) cos t dt). 

The profile will be closed and convex only if a set of N linear 

constraints, r^ 0 for 1 £ i <. N, are satisfied. In practice, CPI 
finds that profile £ which provides the least-squares estimator 
for £ by finding that vector, £, of Fourier coefficients which 
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satisfies the constraints and is as close as possible to the 
vector, of the unconstrained Fourier coefficients Cjv 


v„ (<£) 
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FIGURE 1 


Figure 2 reproduces figures from Ostro and Connelly (1984) 

(i) il lustrating CPI 's geometric conventions and (ii) showing results 
for an asteroid with a D-shaped mean cross section. 
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Fig. 1. Geometry for two-dimensional asteroid light- 
curve inversion. The asteroid is a convex profile rotat- 
ing clockwise, and is shown at rotational phase 6 = 0°. 
The solar phase angle is indicated, as are the aster- 
oid^ illuminated (solid) and unilluminated (dotted) 
portions. The asteroid’s brightness is proportional to 
the orthogonal projection, in the direction of the 
Earth, of the visible, illuminated length 1(8). In the 
normal-angle parameterization of the profile, the point 
P(8) is on the receding (left-hand) limb and the out- 
ward normal n at P{6) points in the positive x direction 
when th^rotational phase is 8. 



Original D-shaped Results of Convex- Profile Inversion 



Rotational Phase 



Fig. 3. Convex-profile inversion of the lightcurve obtained at <b - 30° for a D-shaped, two-dimen- 
sional asteroid, D. Lightcurves, radius-of-curvature functions, and convex profiles are shown in the 
top, middle, and bottom rows, respectively. The left-hand column shows these quantities for D. The 
other columns show results of convex-profile inversion of D's lightcurve when 4> - 30° is assumed 
(middle column) and when <f> » -30° is assumed (right-hand column). The product profile (f) obtained 
using <f> = 30° differs from D because the Fourier representations for the lightcurve and radius-of- 
curvature function (both of which contain discontinuities) have been truncated to 10 harmonics. (See 
end of Section II.C.) The severe distortion in the profile (i) obtained using <£ = - 30° occurs because the 
sign of the solar phase angle is wrong. (See Section II. D.) D’s lightcurve is reproduced as the hand- 
drawn, dashed curve in (d) and (g) for comparison with the model lightcurves (y, dots; computer- 
drawn curves) obtained from the unconstrained and constrained Fourier series x and i respectively. 
In (e) and (h), dots represent r and curves represent k The geometric conventions presented in Fig. 1 
are maintained here. 


FIGURE 2. Two figures from Ostro and Connelly (1984). 
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11 • SUMMARY OF RESEARCH PROGRESS 

The overall objective of the current research is to calibrate and 
apply CPI. Four tasks were specified in the work plan: 

1. Invert high-quality light curves. 

2. Calibrate CPI's sensitivity to departure from ideal geometric 
conditions and to light curve errors. 

3. Explore CPI's sensitivity to non-geometric scattering. 

4. Improve CPI's computational efficiency. 

As the following summaries demonstrate, progress has been 
unexpectedly fruitful in elucidating the power and limitations of CPI, 

f 

in formulating tools for handling profiles as mathematical entities, 
in suggesting novel uses of the profiles, and in developing new 
approaches to testing empirical hypotheses about asteroid shapes and 
scattering properties. 

Inversion of Published .Light-curves 

The entire asteroid light curve literature has been searched and 
catalogued, and lightcurves for some 50 asteroids have been digitized, 
inverted, and subjected to various analyses described herein. Because 
of Eqn. 1, CPI's stability with respect to odd harmonics in r(e) 
deteriorates as § + 0 and -»> 0 for odd n (Fig. 1). A safe rule of 
thumb seems to be to suppress odd harmonics (hence symmetrizing the 
profile) when |(}>| < 10°. The resulting symmetrized profile Pg is 
the best available estimator for the asteroid's symmetrized mean cross 
section Cg. Figure 3 shows CPI output profiles (P or Pg) and 
lightcurves (y, y, y) for seven objects. When I <)>| > 10°, the chosen 
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FIGURE 3. Convex profiles for seven asteroids. In the lightcurve plots, 
large, rectangular symbols are the data y, the solid curve is 
the unconstrained model and the dottea curve is the constrained 
model y corresponding td~the profile. Sources of the lightcurves 
are: Pallas, Binzel (1984) in prep.; Eleonora, Zappala and 

van Houten-Groenveld (1979), A. A. Supp. 35 , 220, Fig. 19; 
Mandeville, Zappala et al . (1983), Icarus 56 , 325, Fig. 28; 
Ra-Shalom, Ostro et al . (1984), Icarus 60 , in press; Ausonia, 
Lagerkvist (1981), A. A. Supp. 44 , 401, Fig. 3; Roxane, Lagerkvist 
et al . (1982), Sun and Planetary System (Fricke and Teleki, eds.), 
p. 209; Hektor, Dunlap et al . (1969), Astron. J. 74 , 796, Fig. 8. 
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sign of <t> corresponds to direct rotation and a southerly view of the 
asteroid. 

Computational Efficiency 

The most time-consuming operation in CPI involves finding that 

A 

Fourier vector £ for the light curve which satisfies the constraints, 
r i 2 l i = 1 to N, and is closest to the unconstrained vector £. 

Our solution of this "quadratic programming problem" uses a 
straightforward, but somewhat inefficient, recursive projection 
al gorithm. 

Lengthy experimentation with various modifications to our 
algorithm suggest that the most economically sound approach is simply 
to pre-order the N constraints from the most negative r^ to the most 
positive r-j . This yields "con vexifi cation" times between 20% and 
50% of those experienced without the heuristic. 

Lightcurves of Geometrically-Scattering Ellipsoids (GSE‘s) 

In efforts to test CPI by inverting lightcurves of known shape 
and scattering law, we realized that whereas it was known that 
lightcurves of GSE's could be generated analytically, a simple formula 
for such lightcurves was unavailable. We derived such a formula and 
discovered properties of GSE lightcurves that let us test the 
hypothesis that a given asteroid lightcurve could be due to a GSE. 
Principal results are: 

1. Given the equation (x T Qx = 1) for an axbxc GSE and 

directions § and § of the Sun and Earth, the lightcurve has 
the form: 
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I = itabc[(^ T Qg) 1 / 2 (s%) 1 / 2 + e T Qs]/2(s%) 1 /2 

2. At opposition (g = §), I 2 has harmonics 0 and 2 only. 

3. The fractional error incurred using an opposition lightcurve 
I 0 to approximate a non-opposition lightcurve I satisfies: 

(I 0 - I)/I £ e s p 2 tan 2 (<t>/2 ) 

where <b is the solar phase angle and p = a/c is the ratio of 
the longest axis to the shortest. Figure 4 plots e vs. <(> for 
several values of p. For p £ 4, the error is less than 1% if 
$ is only a few degrees, but is of order 10% for $ ~ 20°. 

4. The contribution to I of harmonics other than 0 and 2 depends 
on p and <|>. For any lightcurve, we can calculate a statistic 
x , and then use x to find the minimum value of p that permits 
the hypothesis that "the lightcurve can be due to a GSE" to 
stand. The minimum value of p is: 

Pmin * [(1 + x 1 / 2 ) 1 / 2 . l] 1/2 /tan(<D/2) 


where 


x 


^2 2 

l (a + bj 

n*0,2 n " 


2^22 
2a o + Jj (a n + b n> 


< x = (2e + e 2 )? 


where a n and b n are the Fourier cosine and sine coefficients in 
the M-harmonic series for I 2 . Figure 5 plots curves of p m -j n vs. x 
for several values of 41. Positions of asteroid numbers denote minimum 
values of a/c required to forestall rejection of the hypothesis that 
the asteroid is a GSE. 
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Sensitivity of CPI to Viewing/Illumination Geometry 

How reliable is P_ as an estimator for £ when Condition C 
(astrocentric declinations 6 e and 6s of the Earth and the Sun are 
zero) is violated? For any given object, the answer to this question 
depends on 6e, 6s, and <t>. To explore this particular type of 
systematic effect, we have generated and inverted lightcurves of 
GSE's for non -ideal geometries and various values of b/a and c/a which 
sample the parameter space: 10° £ <i> <. 80°, b/a = 0.2 and 0.4, and 

0.05 < c/a < b/a. 

For any given ellipsoid and solar phase angle, the field of 
possible vi ewing/i 11 umi nation geometries is conveniently displayed in 
a plot (Fig. 6) of 6s vs. (6 e - 6s). Figure 6 shows three of 
our several dozen "geometric distortion" plots, selected to illustrate 
a few 2-D cuts through the 5-D space. A major conclusion is that if 
|6 E |, 1 6$| > j<t>j, then distortion is minimal when 5e is between 
|(|>j and -6$ (i.e., when e and s are on opposite sides of the 
equatorial plane) and is most severe when 6s is between |<t>| and 6 e 
( i.e., when s is between e and the equatorial plane). 


Quantitative Comparison of Profiles 

Recognizing the desirability of quantifying the relationship 
between profiles rather than merely making verbal statements of 
subjective impressions, the principal investigator surveyed approaches 
to this type of problem within the pattern-recognition literature. It 
seems appropriate to define the distance Q between any two profiles Pj 
and P 2 as 
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b/a = 0.4 c/a = 0.2 $>=20® b/a =0.2 c/a =0.2 <#>=20° 



S E“ S S 


b/a *0.4 c/a *0.2 $.*40® 



FIGURE 6. "Geometric distortion plots." Application of CPI to GSE lightcurves 
obtained when Condition C ( = = 0) is violated yields profiles 

that deviate from the ellipsoid's mean cross section (£, at the origin 
of each array of profiles). The three arrays illustrate distortion 
for three sets of values of b/a, c/a, and 0. The bottom right figure 
shows some contours of constant $rg or constant . Note that the 

Earth is north of the equator except within the stippled area. 
Distortion is minimal within the region "MIN." 
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Q = 10 min |di/d 10 - d 2 e J0rot /d 2O | 

0 rot 1 

where d-j is the vector of complex Fourier coefficients in the 
expansion of P-j's radius-of -curvature function and d-j 0 is the 
constant, or zero-harmonic term, in that expansion. This distance 
measure is rotation- and scale-invariant and defines a Euclidean 
metric. 

Measures of Noncircularity and Elongation 

We define (i ) the "noncircularity" Q c as the distance between a 
given profile and a circle, and the "maximum breadth ratio" p as the 
ratio of a profile's maximum breadth to its minimum breadth. These 
statistics incorporate all the information in a parent lightcurve, and 
are now being calculated for all our asteroid profiles. 

In the past, the lightcurve amplitude Am has been used as a 
measure of an asteroid's noncircularity and elongation. This approach 
would be justified if every asteroid were a geometrically scattering 
ellipsoid and every lightcurve were obtained at zero solar phase angle 
and under equatorial geometry. For the purpose of using actual 
lightcurves to compare the shapes of real asteroids, the statistics 
Q c and p are a priori preferable to Am and also permit a distinction 
between nonci rcularity and elongation. As a case in point, compare 
the pairs of values of Am, Q c , and p for two lightcurves of 
15 Eunomia at different solar phase angles (Fig. 7), and note values 
of Q c and p for the profiles in Fig. 3, 
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Fig. 7. Convex profiles for asteroid 15 Eunomia. Lightcurves obtained at solar phase angles <j> ** 

-12° (Fig. 8 of Groeneveld and Kuiper, 1954) and <f> = -21° (Fig. 14 of van Houten-Groeneveld and 
van Houten, 1958) were digitized and inverted as described in the text. The signs for <f> correspond to 
retrograde rotation, as deduced by the cited authors. (See Section Il.D.) In each of the top figures, the 
large rectangular symbols are the lightcurve data y , the solid curve represents the lightcurve y derived 
from the unconstrained Fourier fit and the tiny dots represent the lightcurve j derived from the 
constrained fit, i.e., from the convex-profile inversion. In the middle figures, the radius-of-curvature 
functions r and i are represented by symbols and solid curves, respectively. The profiles are shown at 
rotational phase 0 = 0°, and the geometrical conventions given in Fig. I are maintained here. Hence, at 
secondary lightcurve maximum (0 — 320°), each profile would be rotated clockwise —320° and the 
uppermost side would face toward Earth (i.e., toward the top of the page). The relation between the 
rotational phases of the two profiles (and of the corresponding lightcurves and radius-of-curvature 
functions) was chosen to maximize the cross-correlation between 1(0) for = -12° and H0> for <t> = 

-21°. (This approach might prove useful in refining estimates of asteroid rotation periods.) In these 
and all other inversions presented in this paper, each Fourier series was truncated to M m 10 harmon- 
ics and the inequality constraints (8) were enforced at N =* % equally spaced values of 0. 

FIGURE 7. Lightcurve amplitude Am-, noncircularity fl c , and maximum breadth 
ratio for two independent lightcurves of 15 Eunomia have been 
added to this reproduction of Fig. 7 from Ostro and Connelly (1984). 
Note also the average of the two profiles, and its values of jl and ft. 

C 
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Algebraic Manipulation of Profiles 

To the extent that ideal conditions are satisfied, profiles 
obtained for different solar phase angles should be similar. In 
different words, they should not be "far apart" because each is an 
estimator for the asteroid's (unique) mean cross section. Therefore 
it makes sense to obtain appropriately weighted averages of 
independent profiles. In practice, this is most easily done by 
manipulating the Fourier vectors £jj for r-j ( 9 ) corresponding to the 
Pj . Thus, CPI is able to combine the information contained in any 
number of independent lightcurves to yield a 2-D average of an 
asteroid's 3-D shape. Figure 7 shows the average of two Eunomia profiles. 

Sensitivity of CPI to Non-Geometric Scattering and Non-Convexity 

Efforts are underway to calibrate the distortion introduced into 
a profile P_ as an estimator for the mean cross section £ when ideal 
conditions A (convexity) and B (geometric scattering) are violated. 

J. Lambert (Univ. of New Mexico) has used a numerical integration 
routine to generate lightcurves of ellipsoids having non-geometric 
scattering laws, for various geometries. His lightcurves, which 
include a few "standard" GSE's, are being inverted at Cornell. First 
results indicate that limb darkening results in a profile that is more 
elongated than the parent ellipsoid's mean cross section. In other 
words, (Ms " less circular" than C_. In contrast, departure from ideal 
viewing/i 11 umi nation geometry (Condition C) frequently causes £ to be 
more circular than C (Fig. 6). Perhaps these two systematic effects 
offset each other. In any event, given adequate prior knowledge of 
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viewing/illumi nation geometry or scattering law, we can try to correct 
P_ for non-ideal conditions via ad hoc addition or subtraction of some 
fraction of a circle. 

In a different experiment, M. A. Barucci (Instituto de 
Astrofisica Spaziale, Co.ns1.glio Nazionale delle Ricerche, Rome) has 
generated laboratory lightcurves of a stone fragment for a variety of 
viewing/illumi nation geometries. The lightcurves have been inverted 
at Cornell, and photographs of the fragment have been requested from 
Barucci to provide a "ground truth" fiducial for calibration of CPI. 

f " 

Propagation of Lightcurve Noise into Profile Error 

The previous paragraphs have focused on systematic sources of 
error in P» i.e., on distortion in P_ caused by departure from ideal 
conditions. Since actual lightcurves are contaminated with 
measurement errors, any profile obtained via CPI may be corrupted 
accordingly, and it is necessary to examine the nature of this error 
propagation. 

Whereas actual lightcurve measurement errors may arise from 
systematic sources (e.g., sky-background removal, intermittent cirrus, 
imperfect color matches between asteroid and comparison star, etc.) as 
well as from photon-counting stochastic errors, one can always compute 
an "effective" standard deviation in the overall noise due to all 
sources. As this standard deviation becomes an increasingly large 
fraction of the lightcurve amplitude, P^will become increasingly 
distorted. The distortion will be magnified if the lightcurve 
contains big gaps or is sampled at an insufficient rate. 
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We have investigated the three effects of noise, gaps, and 
undersampling, acting alone and in concert, on profiles obtained from 
lightcurves of GSE's and the D-shaped object in Fig. 2. Propagation 
of statistical error is being modeled by using a random number 
generator to produce Gaussian noise samples, scaling the standard 
deviation to some fraction of the lightcurve amplitude, adding the 
noise to the lightcurve, and doing the inversion. This entire process 
is repeated to build up a family of realizations of the random process 
in the form of a set of profiles, Pj. The distortion in the Pj is 
evident visually and can be quantified by the distances of the Pj 
from C. Figure 8 illustrates some results of five-profile simulations 
on GSE lightcurves. (Distortion is worse for the D-shaped object.) 

At this point, an appropriate guideline seems to be that 
lightcurves targeted for CPI should be sampled at intervals no larger 
than 5° and should have noise levels < 3% of the lightcurve 
amplitude. Since CPI's ideal conditions on viewing/illumination 
geometry, etc., will rarely be satisfied exactly, it clearly is vital 
that the lightcurve data be of high quality. 

Current error-propagation work is focusing on actual asteroid 
lightcurves, with the objective of placing formal errors on Q c and p 
as well as on the profiles themselves. 



FIGURE 8. Error propagation into profiles obtained from lightcurves of a 
GSE with b/a - 0.4, c/a = 0.2, and 6 e = 6s - Q.* Sets of 
five simulation profiles are labeled with the sampling inter va l 
in degrees, the noise level (expressed as a percentage of the 
lightcurve amplitude), and the distance Q from the parent 
ellipsoid's mean cross section. "INTER -EX OUT" denotes 
deletion of four, 40° inter-extremum sections of data, and "MIN 
OUT" denotes deletion of a 40° section of data centered on one 
minimum. 
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III. FUNDING STATUS 

Spending during the current report period was at the anticipated 
rate. A graduate student in Applied and Engineering Physics, Mark 
Dorogi, was hired full time over the summer as a research assistant. 

A Physics senior, Anthony Ferro, will be employed as a part-time 
programmer this fall. 

IV. PUBLICATION OF RESULTS 

The following papers report results of the research supported 
under NAGW-543: 

1. Convex Profiles from Asteroid Lightcurves, S. J. Ostro and 

R. Connelly, Icarus 57, 443-463 (1984). 

2. Ellipsoids and Lightcurves, R. Connelly and S. J. Ostro, 

Geometriae Dedicata , in press (1984). 

A paper (abstract appended) reporting some of the more recent 
results has been submitted to the 1984 Meeting of the Division for 
Planetary Sciences of the American Astronomical Society. 

Additionally, the principal investigator has given seminars on 
CPI at the Arecibo Observatory, the Planetary Science Institute 
(Tucson), the Jet Propulsion Laboratory, and Cornell. 
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Abstract of paper submitted to the 16th Annual Meeting of the Division 
for Planetary Sciences of the American Astronomical Society: 


Convex-Profile Inversion of Asteroid Lightcurves : 
Calibration and Application of the Method 

S. J. Ostro, M. D. Dorogi, and R. Connelly (Cornell) 

Convex-profile inversion (Ostro and Connelly* 1984, 
Icarus 57 , 443) is a method for obtaining a convex 
profile, £, from asteroid lightcurve data. Whenever 
four ideal conditions are satisfied, £ is an estimator 
for the asteroid's "mean cross section" £, a convex set 
defined as the average of all cross sections £(z) cut 
by planes a distance z above the asteroid's equatorial 
plane. £ is therefore a 2-D average of the asteroid's 
3-D shape. The ideal conditions are: (A) each curve 

£(z) is convex, (B) the asteroid's scattering law is 
uniform and geometric, (C) the astrocentric 
declinations 65 , 5e of the Sun and Earth are zero, 
and (D) the solar phase angle <j> * 0. The method has 
been tested by inverting lightcurves generated 
analytically for geometrically scattering ellipsoids 
(GSE's) with semiaxes a £ b £ c. Using a suitably 
defined "distance measure" to quantify the difference 
between any two profiles, we have calibrated the 
deviation of £ from £ for GSE’s as a function of 
lightcurve noise level, rotation-phase sampling 
interval A9, and departure from ideal conditions. The 
coupling between these factors is considerable. If the 
rms noise due to all sources is ~ 3% of the 
lightcurve's peak-to-valley amplitude, then £ « £ if 
A9 < 5°. If |6e|» I | > then distortion of 
£ is minimal when 6£ is between 0 and -65 but 
severe when 6s is between 0 and 6 e* We currently 
are calculating £ for high-quality, published 
lightcurves. The "distance" between £ and a circle 
provides a gauge of the asteroid's nonsphericity and 
(in contrast to the lightcurve amplitude) incorporates 
all the information contained in the lightcurve. This 
research was supported by NASA Grant NAGW-543 (SJO, 

MDD) and NSF Grant MCS77-01740 (RC). 



